Non-trivial fixed point in a twofold orbitally degenerate Anderson impurity model 
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We study a twofold orbitally degenerate Anderson impurity model which shows a non-trivial fixed 
point similar to that of the two-impurity Kondo model, but remarkably more robust, as it can only be 
destabilized by orbital- or gauge-symmetry breaking. The impurity model is interesting per se, but 
here our interest is rather in the possibility that it might be representative of a strongly-correlated 
lattice model close to a Mott transition. We argue that this lattice model should unavoidably 
encounter the non-trivial fixed point just before the Mott transition and react to its instability by 
spontaneous generation of an orbital, spin-orbital or superconducting order parameter. 
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PACS numbers: 

When a metal is driven by strong correlations towards 
a Mott insulator, an incoherent component of the single- 
particle spectrum slowly moves away from the quasipar- 
ticle resonance and smoothly transforms into the Mott- 
Hubbard side-bands. Analogous behavior is displayed by 
an Anderson Impurity Model (AIM) from the mixed va- 
lence to the Kondo regime. This is suggestive of similar 
physical processes underlying the dynamics of AIM's and 
of strongly-correlated electron systems across the Mott 
metal-to-insulator transition (MIT), even though a rig- 
orous relationship holds only in infinite dimensions, as 
shown by Dynamical Mean Field Theory (DMFT). Fur- 
thermore, DMFT shows that the Mott-Hubbard bands 
split off from the quasiparticle resonance quite before the 
MIT occurs lj, suggesting that it is rather the physics 
of the AIM in the Kondo regime to be significant near 
the MIT. In that limit, the AIM maps onto a Kondo 
model with the number of conduction channels always 
such as to perfectly screen the impurity. For that reason 
one would exclude that the appealing non-Fermi liquid 
physics of the overscreened multi-channel Kondo effect 
may ever appear close to the MIT. 

That expectation is partly wrong. In this Letter we 
study the phase diagram of an AIM which does contain 
a non-trivial fixed point. We show that the physical be- 
havior around this fixed point resembles that displayed 
by the two-impurity Kondo model. We also argue that 
the lattice model, where the physics of the above AIM 
should be relevant, necessarily encounters this fixed point 
just before the MIT, and discuss possible consequences. 

We consider the two-orbital AIM Hamiltonian: 

HAIM = [ £fc C l,aa C k,aa + ( V k4,aa d aa + H.C. 

k.a.cr 
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where c\ are conduction-band creation operators in 



orbital a = 1, 2 and spin a, d\ a impurity ones and = 
E aCT <4r<C- For H K = 0, the AIM has a large SU(4) 
symmetry which is lowered down to a spin SU(2) times 
an orbital 0(2) by a Hund's rule- like coupling: 
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where T* = ± £ CT 
are orbital pseudo-spin operators with r 
matrices. For later convenience, we 
duce the impurity spin, S l — \ £ a £ Q 



= x,y,z, 
s the Pauli 
also intro- 



a aa 



a/3 d a /3, 



and the impurity spin-orbital operators W v 

> J2a,b Ea,(3 d aa <T l a/3 T^ b d b/3 . 

In the Kondo regime, U much larger than the conduc- 
tion bandwidth, two electrons get trapped by the im- 
purity in a configuration identified by total spin, S, to- 
tal pseudo-spin, T, and their z-components, with energy 



E(S,S Z ;T,T Z ) = K 



By Pauli prin- 



T(T+1) - {T Z Y 
ciple, two electron configurations have either S = 1 and 
T = 0, or S = and T = 1. If K > 0, the lowest energy 
configuration has 5 = 1 and T = 0, the conventional 
Hund's rules. The impurity behaves effectively as a spin 
S=l which may be Kondo-screened by the two conduc- 
tion channels (5 = tt/2 phase shift). On the contrary, if 
K < the non-degenerate S = T = 1 and T z = state 
has lowest energy. Here we do not expect any Kondo 
effect, i.e. S = 0. This situation is analogous to the two 
S=l/2 impurity Kondo model (2IKM) in the presence 
of a direct exchange between the impurities. There, it 
is knowndHQ that under particular circumstances an 
unstable fixed point (UFP) separates the Kondo-screened 
regime from the one in which the two-impurities couple 
together into a singlet. In our model (JIJ, that circum- 
stance is realized thanks to the 0(2) orbital symmetry, 
as shown later, hence an analogous UFP should exist. 

We have studied the AIM (JJJ in the Kondo regime 
by Wilson's Numerical Renormalization Group (NRG), 
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FIG. 1: NRG flow of the lowest energies of the levels labeled 
by {Q, T z ,S), where Q is half of the added charge. The right 
and left flows correspond to a relative deviation SK,/K* — 
±4 • 10~ 3 from the fixed point value K*, respectively. The 
phase diagram is sketched in the inset. 



closely following the original work by Jones and Varma 
for the 2IKM[J0|. We have also developed a comple- 
mentary analysis based on abelian bosonization, which 
provides a better characterization of the UFP. 

In the inset of Fig. 1 we sketch the phase diagram of 
||TJ in the Kondo regime, as determined by the flow of the 
low energy spectrum obtained by NRG. At fixed Kondo 
exchange we find indeed an UFP upon varying A. For 
K > A*, the model asymptotically flows to a Kondo- 
screened fixed point, see right panel in Fig. 1, while for 
A < A", < it flows towards a non-Kondo-screened fixed 
point, see left panel. The intermediate crossover region, 
also shown in Fig. 1, identifies the UFP [3. We notice 
that CJ has a larger impurity Hilbert space then the 
2IKM, which contains, besides the 5=1 and the T z = 
5 = configurations, also the T z = ±1 5 = dou- 
blet, absent in the 2IKM. In spite of that, the low-energy 
spectra at the UFP's are the same for both models. In 
Fig. 2 we plot the ground state average values of the im- 
purity operators 5 2 , T 2 and (T z ) 2 . For large and positive 
A, the impurity freezes into the 5 = 1 T = configu- 
ration while, for very negative values, into the T z = 
5 = one. At the UFP (5 2 ) = 1/2, (T 2 ) = 3/2 and 
((T z ) ) = 1/4. To prove that our UFP is connected with 

the 2IKM one, we have added to (JTJ a term G z (T Z J , 

with G z > 0, which pushes upward the energy of the 
T z = ±1 5 = doublet absent in the 2IKM. In the in- 
set of Fig. 2 we plot (5 2 ) as function of G z at the UFP, 
whose position depends on G z too. We do find that (5 2 ) 
smoothly reaches the 2IKM unit value for large G z . 

The approach to the two stable fixed points, A < A* 
and K > A* , can be described by the local perturbation 
left behind by the impurity which has either disappeared, 




FIG. 2: Average impurity quantum numbers as function of 
K. In the left inset the behavior around the UFP is shown 
in more detail, while in the right inset (S 2 ) along the path 
towards the 2IKM is displayed (see text). 



S = 0, or been absorbed by the conduction sea, 6 = it/2: 
H. = -Y / t*(fl a J haa + H.c.)+^(n -2) 2 



(3) 



where "0" labels the first available site of the Wilson 
chain, i.e. the actual first site for A < A*, 6 = 0, and 
the second site for A > A*, 6 = tt/2. Numerically we 
find Jt* ~ 2Z7* ~ 2Js* ~ 32i* — * oo upon approaching 
the UFP, see Fig. 3. The diverging t» implies a singu- 
lar (oc (A — A*)~ 2 ) impurity-contribution to the specific 
heat, just like in the 2IKM. Additional information are 
provided by the Wilson ratios R4 = (c v Sxi) / (Xi$ c v), 
where i = S,T Z refer to spin and orbital susceptibili- 
ties, i.e. to those response functions related to conserved 
quantities, hence accessible by Fermi liquid theory 0, 0] ■ 
Rs and Rt* are shown in Fig. 4 and appear to vanish 
at the UFP. By analogy with the 2IKM, there are two 
other susceptibilities which arc instead expected to be 
singular: the susceptibility Xst* to a field which couples 
to the relative spin operator W lz , i = x,y,z, (the stag- 
gered spin-susceptibility in the 2IKM), and the pairing 




FIG. 3: ^-dependence of the effective couplings in Eq. (J^J. 
In the inset the region around K = is shown. 
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FIG. 4: Spin (Rs) and orbital (Rt*) Wilson ratios as func- 
tions of K. Notice that the K — 0, SU(4)-point, as well as 
large K, S — 1 impurity (shown in the inset), values coincide 
with known results. 



susceptibility xsc m the Cooper channel c^Cg. + cLcj, . 
Those are not accessible by Fermi liquid theory. Yet one 
can get a rough estimate of them by the corresponding 
scattering amplitudes at zero external frequencies. They 
are given respectively by Vst* = — 2E/» + Js* — Jt* and 
Tsc — 2C/* — 3 Js* — Jt*, hence are negative (correspond- 
ing to an enhancement of the response) and diverge sim- 
ilarly approaching the UFP. The physics underneath is 
the same of the 2IKM, and has been exhaustively dis- 
cussed in Rcf. 4J. The UFP has a residual entropy 
ln-\/2. Away from the UFP, this entropy is quenched 
below a temperature scale T* ~ \K — K*\ 2 , implying a 
specific heat coefficient 7 ~ 1/T*. The rest of the impu- 
rity entropy is quenched at higher temperatures of order 
T K ~ \K\. At the UFP, 7 is finite, whilst both xst* and 
Xsc display a | lnT| singularity. 

The stability of the UFP is more easily accessed by 
abelian bosonization, following Ref. [8j on the 2IKM. In 
the large U limit, maps onto the Kondo model 

3 
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where a 1 , f % and uj ij are, respectively, the conduction- 
electron spin, orbital and spin-orbital densities at the 
impurity site. As usual in abelian bosonization we al- 
low for anisotropy: Jf = Jj 7^ JJ, and similarly for J T 

Jw> for hJ ^ z - Tne 



and Jyy- We further assume Jyy 



spin ' 



anisotropic Kondo model 10} has a continuous 0(2) 
O(2) orbita ;x \J(l) ch arge symmetry, which is useful to de- 
compose into U(l) spin xU(l) or bj to ;XU(l) fc arge plus two 
discrete symmetries: (i) a IIJ pm rotation: c a (T <-» c a ,-o 
and d a ,a <-* d a ,~a', and (ii) a H° rb rotation: o\ a <-> 
C2a and dxa *-* d^. By abelian bosonization|9j, we 
write the s-wave scattering components of conduction 



electrons as chiral one-dimensional fermions c a . a (x) — 
1/v2tt a F Qj(T exp [— i<f> ajCr (x)\, where 4>a,a are chiral free 
Bose fields, a a short distance cut-off, and the Klein fac- 
tors -F a ,(T are Grassman variables enforcing proper anti- 
commutation relations. Next, we introduce the combina- 
tions: cj) c = (</>i T + cf>n + (j> 2 ] + 02j,)/2, <j) s = (0i T - cj>n + 
^2T-02i)/2, 4>f = (0it — 02t — 02i)/2, (/) s f = (0it~ 
0U ~~ 02T + 02i)/2. After applying the canonical transfor- 



mation exp 



exp 



if z 4>f(0) , Eq. |@J| can be re- 



fermionized via ^(x) = \j \J1-kol Ft exp [— i(j)b(x)], where 
b = c, s, /, s/|9j. For a particular value of JJ = Jf., the 
end result is an effective model where only ty s f is coupled 
to the impurity, just like in the 2IKM0. 

To locate the UFP, we follow the same strategy of 
Ref. : we assume K large compared to the conduction- 
bandwidth and search for an accidental ground state de- 
generacy in that part of the effective Hamiltonian involv- 
ing just the impurity and the i^'s: 



imp 



= H K + X S (S Z ) 2 + X T (f zXZ 



7-L 

2na 



F$Ft W— + F f F^ W-+ + H.c] .(5) 



As and At are cut-off dependent functions of J s , J T . 
For a specific K* < 0, we find that the impurity state 
I 0) = I S = 0, S z = 0; T = 1, T z = 0) is degenerate with: 

\l) = S^.(^ Ff |0,0;1,+1)+F; I 0,0; 1,-1)) 

+ ^(F S |1,+1;0,0)-F] |1,-1;0,0)), (6) 

where 9 depends on the Hamiltonian parameters. For 
our model QJ, 9 should be equal to 7r/4 to reproduce 
the observed UFP average values of S 2 , T 2 , and (T z ) 2 . 

If we added the term G, (f z ) , 9 should increase with 

G z , reaching the 2IKM value of 9 = 7r/2 for large G z . 
The Klein factors in © show that | 0) and | 1) differs 
by one fermion, justifying the introduction of a fictitious 
fermion connecting that doublet: f'\0) = \l). 

The low-energy Hamiltonian close to the UFP, Hjjpp, 
is then obtained by projection onto the above doublct- 
subspace. Including up to dimension 3/2 operators, 



Hufp — Ho + Ao 



*L(0)-* sf (0)l (p + f ) (7) 



Ai^ 



*^(0)-* fl/ (0)l (p-f)+SK*pf 



with H the free Hamiltonian for the fj^j's, and 5K* 
the deviation from the fixed-point value K*. Xq and Ai 
are model dependent parameters. As expected, Eq. Q 
has the same form as in the 2IKM@. The UFP Hamil- 
tonian [first line of Eq. |Q] is a resonant level model in- 
volving one Majorana fermion * s f — ^sf hybridising with 
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ft + f. The combination ft — f is free and is respon- 
sible for the ln\/2 UFP residual entropy. The relevant 
term (dimension 1/2) proportional to 5K* describes the 
deviation from the UFP, while the Ai-term is the leading 
irrelevant operator (dimension 3/2). Other possibly rele- 
vant operators are instead not allowed by the symmetry 
properties of |@J, which have to be preserved by Hufp 
too. For instance, among the particle-hole symmetry 
breaking terms allowed in the2IKM[|, only the marginal 
one, which does not spoil the UFP properties, may ap- 
pear in our model, since the relevant operator, bosoniza- 

tion of which is given by 1^(0) + * / (0) (/ - ft)^, is 

here forbidden by ~U(l)orMtai- In fact, while / is invari- 
ant under a V(l) or bitai rotation parametrized by a phase 
a, due to the Klein factors in ©, 9f transforms into 
e 2la ^>f. Indeed all relevant perturbations which destabi- 
lize the UFP correspond to physical instabilities of model 
||TJ, unlike what happens in the 2IKM. For instance, 
the relevant terms [\P C (0) ± *J(0)] (/ - ft), of dimen- 
sion 1/2, break \5{l) c harge- Therefore, gauge symmetry 
breaking destabilizes the UFP, which explains the singu- 
lar behavior of xsc- Analogously, xst* is the response 
to a field which breaks SU(2) sp i n x H° rb and allows the 

relevant termsjl^ 1*1/(0) + * a/ (0)l (/ - / f )|l @ and 

± * s (0)] (/ - ft) . Besides those two susceptibili- 
ties, also XT a and XST" , with a — x,y, are logarithmically 
diverging, being related to fields breaking XJ(l) or bital- 

We now turn to our original motivation and discuss the 
possible relevance of the above results to the physics of 
the Mott transition. Take a lattice model with an on-site 
interaction of the same form as in CJ-©, with inverted 
Hund's coupling K < 0. This may occur if the electrons 
are Jahn- Teller coupled to two degenerate weakly disper- 
sive optical phonons by g ^2r 

(qi R fl + q 2R f£), where 
qm are the phonon coordinates on site R. This cou- 
pling gives rise to a retarded electron-electron interaction 
which reduces to J5J with K ~ —g 2 /ojo when the typical 
phonon frequency luq is much larger than the quasiparti- 
cle bandwidth. Alternatively, two single-band Hubbard 
planes/chains coupled by JJ^rSir. ■ S2R, where 1 and 
2 refer to the two planes/chains and J > 0, would also 
display a similar behavior. 

When K = the lattice model should undergo a MIT 
at some finite U c in the absence of nesting. If \K\ <C U c , 
the physics of the metallic phase near the MIT should 
resemble that of the AIM, Eq. QJ, in the Kondo regime. 
Since the width of the quasiparticle resonance, i.e. the 
effective Kondo temperature Tk, vanishes at the MIT, 
Tk ~ U c — J/Q, the system is forced to enter the crit- 
ical region around the unstable fixed point, \K\ ~ Tk, 
before the MIT occurs. However the instability of the 
AIM around the UFP towards the orbital 0(2) or charge 
U(l) symmetry breakings should transform in the lat- 
tice model into a true bulk instability. Namely, at least 



within DMFT, we expect that the self-consistency condi- 
tion which relates the impurity Green's function with the 
local Green's function of the bath enlarges the UFP into 
a whole region where the model undergoes a spontaneous 
symmetry breaking. This would open up new screening 
channels for those degrees of freedom which survive be- 
low T K down to T* - \T K + K\ 2 /T K (K < 0) and are 
responsible of the finite entropy at the UFP. If the band- 
structure lacks nesting or Van Hove singularities, orbital 
or spin-orbital instabilities are not competitive with the 
Cooper instability 0. This suggests a superconducting 
region just before the MIT, which would be remarkable 
since the bare scattering amplitude in the Cooper chan- 
nel is U + K, hence repulsive for U ~ U c >■ \K\. We 
believe that this phase is analogous to the strongly cor- 
related superconductivity recently identified by DMFT in 
a model for tetravalent alkali doped fullerenes^2$- The 
latter model maps by DMFT onto a threefold degenerate 
AIM with inverted Hund's rules, mimicking a t (g> i? dy- 
namical Jahn- Teller effect. Although different from our 
model 01-101, it contains the essential physics we have 
described in this work; namely the competition between 
the Kondo- and an intra-impurity-screening mechanism. 
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Under U° x r 
* c <-> -vjr c 

The Bethe-Salpeter equations for bulk susceptibilities de- 
pend not only upon the irreducible vertices, which are 
local in DMFT, but also on products of lattice Green's 
functions in the appropriate particle-particle (p-p) or 
particle-hole (p-h) channels. While this product in the p- 
p channel is always singular by the Cooper phenomenon, 
in the p-h channels singularities may arise only due to 
band structure effects. Therefore, even though local sus- 
ceptibilities hence irreducible vertices are singular in p- 
p and p-h channels, the bulk p-p susceptibility is more 
likely the dominant one. 
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